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Over  two-thirds  energy  of  fuel  consumed  by  an  automobile  is  discharged  to  the  surroundings  as  waste 
heat.  The  fuel  usage  can  be  more  efficient  if  thermoelectric  generators  (TEG)  are  used  to  convert  heat 
energy  into  electricity.  In  this  study,  a  thermoelectric  module  composed  of  thermoelectric  generators  and 
a  cooling  system  is  developed  to  improve  the  efficiency  of  an  IC  engine.  Two  potential  positions  on  an 
automobile  are  chosen  to  apply  this  module,  e.g.  exhaust  pipe  and  radiator  to  examine  the  feasibility.  To 
predict  the  behaviors  of  this  module,  a  one  dimensional  thermal  resistance  model  is  also  build,  and  the 
results  are  verified  with  experiments. 

The  maximum  power  produced  from  the  module  is  51.13  mWcm'2  at  290  °C  temperature  difference. 
The  model  results  show  that,  TE  module  presents  better  performance  on  the  exhaust  pipe  than  on  the 
radiator. 

©  2009  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

According  to  the  Otto  cycle,  during  the  operating  process,  it  is 
necessary  for  internal  combustion  engine  to  discharge  heat  to  the 
atmosphere  for  completing  the  thermodynamic  process.  Conse¬ 
quently,  only  about  30%  energy  released  from  consumed  fuel  is 
converted  as  propelling  force,  the  other  70%  is  either  discharged  by 
exhaust  gas,  or  expelled  by  the  cooling  system.  Energy  lost  from 
dispersed  heat  can  be  recovered  by  several  ways,  thermoelectric 
generator  (TEG)  is  considered  useful  by  several  auto  makers. 
Through  the  effect  of  temperature  difference,  waste  heat  can  be 
converted  to  electricity  from  the  exhaust  pipe  and  radiator  on  an 
automobile. 

In  1988,  Birholz  et  al.  [1]  presented  the  first  TEG’s  application  on 
the  automobile.  In  this  research,  a  single  TEG  unit  using  FeS^  as 
material  was  adopted  to  produce  1  W  electric  power.  Bass  et  al. 
[2-6]  applied  an  array  with  72  pieces  of  TEG  on  a  diesel  truck.  By 
maintaining  230  °C  and  30  °C  at  hot  and  cold  sides  of  TEGs 
respectively,  energy  conversion  efficiency  of  4.5%  was  achieved. 
Kobayashi  et  al.  [7]  used  the  same  amount  of  SiGe  TEGs  on 
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a  3000  c.c.  petrol  engine  vehicle.  A  rectangular  section  of  exhaust 
pipe  was  used  as  the  hot  source  and  attached  by  one  side  of  TEG.  On 
the  other  side  of  TEG,  a  water  cooling  system  was  used  as  a  heat 
sink.  At  60  km/h  vehicle  speed,  the  exhaust  temperature  was 
1141  °C,  and  the  temperature  difference  on  TEG  was  563  °C.  The 
maximum  power  generated  from  a  single  TEG  was  1.2  W  at  0.7  V 
with  0.9%  thermal  efficiency.  Ikoma  et  al.  [8]  applied  HZ-14s 
(Bi2Te3)  on  a  petrol  engine  vehicle.  The  maximum  power  output 
from  TEGs,  at  vehicle  speed  of  60  km/h,  was  193  W  with  2.9% 
conversion  efficiency.  Thacher  et  al.  [9]  attached  HZ-20s  (B^Tes)  to 
a  light  truck  engine.  The  experimental  results  showed  the  power 
output  from  TEGs  increased  with  engine  speed,  and  the  whole 
system  can  produce  up  to  330  W  power. 

Matsubara  [10,11]  used  two  materials,  YbFe3CoSbi2  and 
YbFe3.6Ni0.4Sbi2,  to  develop  a  new  TEG.  16  of  them  were  stacked  on 
the  exhaust  pipe  of  a  2.0  L  passenger  car.  A  maximum  power  output 
of  266  W  was  gained  under  475  °C  temperature  difference  condi¬ 
tion.  Recently,  Champier  and  co-workers  [12]  applied  TE  modules 
on  biomass  cook  stoves  to  generate  electricity  to  power  the  fan  and 
light.  This  commercial  TEG  is  made  of  Bismuth  Telluride,  the  same 
material  used  in  this  work. 

Many  experiments  have  been  done  to  apply  TEGs  on  an  auto¬ 
mobile,  however,  only  a  few  works  simulated  the  behaviors  of  TEG 
and  TE  module.  The  presented  model  helps  to  understand  the 
characteristics  of  TEG,  and  the  effects  of  engine  speed  and  coolant 
temperature  of  radiator  on  the  TE  module. 
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Heat  Source 


Fig.  1.  Schematic  of  a  working  thermoelectric  generator. 


2.  Experimental 

2.1.  Basics  of  TEG 

The  working  principles  of  TEG  include  Seebeck  effect,  conduc¬ 
tion  effect  and  Joule  effect.  Seebeck  effect  explains  the  phenom¬ 
enon  that  electric  current  induced  by  temperature  difference  in  an 
electrical  conductor,  which  simply  describes  the  main  function  of 
a  TEG.  Typically,  a  commercial  TEG  is  composed  of  ceramic 
substrates,  electrical  insulators,  electrical  conductors,  and  N-type/ 
P-type  semiconductor  block.  By  applying  different  temperature  on 
both  sides  of  it,  electric  current  is  induced  because  of  the  thermo¬ 
electric  effect.  Fig.  1  shows  the  schematic  of  a  working  TEG. 

In  Fig.  1,  the  rates  of  supply  heat  (Qh)  and  removal  heat  (Qc)  can 
be  calculated  at  the  hot  and  cold  junction  respectively. 

Qh  =  aITH  +  K(TH-TL)-ll2Rc  (1) 


the  upper  heating  surface.  10  Thermocouples  were  installed  on  the 
block  to  measure  the  temperatures.  By  collecting  the  data  of 
temperature  differences  between  the  inside  and  outside  positions 
of  block,  the  heat  loss  can  be  determined  by  Fourier’s  Law  of  Heat 
Conduction.  The  total  heat  loss  (Qi)  consists  of  four  side  losses  (Qs) 
and  one  base  loss  (Qb). 

Qi  =  Qs  +  Qb  (3) 

Qs  =  4Ms  T,  s  ~  To  s  (4) 

LS 

0b  =  Mb  T',b  7  —  (5) 

Lb 

2.3.  Performance  coefficients  of  TEG 

The  performance  of  TEG  can  be  characterized  with  Seebeck 
coefficient  (a:),  internal  electrical  resistance  (Rq)  and  thermal 
conductance  (K). 

The  power  generated  from  a  TEG  can  be  measured  either  by 
multiplying  its  output  current  with  voltage,  or  observing  the  vari¬ 
ation  of  heat  flux. 

Po  =  V0/0  =  Qh  -  Qc  (6) 

Subtract  equation  (2)  from  equation  (1) 

Po  =  «/(Th  -  Tl)  -  I2Rc  (7) 

the  voltage  can  be  expressed  as 

Vo  =  ^  =  a(TH-TL)-IRc  (8) 

1 0 

when  the  power  generation  occurs  at  the  condition  of  open  circuit 
(1  =  0),  equation  (8)  becomes 

Vo  =  a(TH  -  Tl)  (9) 

The  Seebeck  coefficient  (a)  can  be  derived 


Qc  =  «/Tl  +  K(Th-Tl)+1/2Rg 


(2) 


2.2.  TEG  performance  tests 

Two  different  brands  of  commercial  TEGs,  H-type  and  T-type, 
were  adopted  for  the  performance  tests.  Table  1  shows  their 
geometry  characteristics.  The  heat  source  to  the  hot  side  and  the 
cooling  method  for  the  cold  side  will  influence  the  performance.  A 
specific  temperature  difference  on  TEG  is  required  to  maintain 
a  better  outcome.  Hence,  a  heater  was  used  to  provide  continual 
energy  to  the  hot  side,  and  a  micro  channel  heat  sink  was  attached 
to  the  cold  side  to  take  the  heat  energy  away.  Water,  as  a  coolant, 
circulated  between  the  heat  sink  and  a  low  temperature  circulator 
bath.  Fig.  2  presents  the  structure  of  experimental  setup. 

The  heater  was  placed  into  a  hollow  Teflon  block  to  prevent  the 
heat  loss  from  the  ambient.  The  block  enclosed  the  heater  except 


_  Vo_ 

(Th  ~  7l) 


(10) 


Apply  equation  (10)  to  equation  (8)  to  obtain  the  value  of  internal 
electrical  resistance  (Re).  Consequently,  substitute  a  and  RG  into 
equation  (1)  to  derive  the  value  of  thermal  conductance  (I<).  Take 
the  first  order  partial  derivatives  of  equation  (7)  with  respect  to  I.  By 
assuming  the  result  on  the  left  hand  side  of  equation  is  zero  [13], 
the  output  current  J0  becomes 


,  _  «(Th  -  TL) 

'°  "  2RG 

The  thermal  efficiency  of  TEG  can  be  calculated  by 


(id 


(12) 


2.4.  Heat  flux  analysis  of  the  TE  module 


Table  1 

Geometry  characteristics  of  TEG. 


Model 

Length  (mm) 

Width  (mm) 

High  (mm) 

Weight  (g) 

HZ-2 

29 

29 

5.08 

13.5 

TGM-127 

40 

40 

3.9 

32 

TEGs  with  a  heat  exchanger  system  were  designed  to  recover  the 
waste  heat  from  two  positions,  e.g.  exhaust  pipe  and  radiator  on  an 
automobile.  To  better  understand  the  behaviors  of  each  component 
in  the  module,  a  mathematic  model  was  built  to  analyze  the  change 
of  heat  flux  and  predict  the  performance  of  system. 
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Fig.  2.  The  integral  system  for  TEG  performance  test. 


2.4 A.  Thermal  resistance  model  of  the  waste  heat  recovered  from 
exhaust  pipe 

Fig.  3  illustrates  a  thermal  resistance  model,  which  recovers  waste 
heat  from  the  exhaust  pipe.  This  model  can  be  classified  into  three 
sections,  the  hot  side,  the  cold  side  and  TEG  itself.  The  surrounding 
fluid  of  the  hot  side  is  exhaust  gas  and  the  cold  side  is  the  ambient  air. 
The  heat  energy  of  exhaust  gas  depends  on  its  temperature,  however, 
the  temperature  increases  with  the  engine  speed. 

Thermal  resistance  of  hot  side  is  composed  of  interface  thermal 
resistance  (R[),  material  thermal  resistance  (Rm)  and  convectional 
thermal  resistance  (KConv)-  The  thermal  resistance  of  cold  side  is 
composed  of  interface  thermal  resistance  (/%),  spreading  thermal 
resistance  (Ksp),  material  thermal  resistance  (Rm)  and  convectional 
thermal  resistance  (KConv)- 

When  two  materials  contact  with  each  other,  thermal  resistance 
occurs  at  the  interface.  It  can  be  expressed  as 


Ri  = 


tj 

Me 


(13) 


There  is  spreading  thermal  resistance  when  heat  flows  from 
a  small  area  to  a  bigger  area  [14] 


Rsp 


7Wb 

Qc 


—  ^sp,2  ^sp,l 


where 


Rsp,l 


\/Xc  x  Q 

kb  'JkXTAg 


AkbAb^conv  +  tanh(Atb) 
1  +  A/<bAb^convtcmh(Atb) 


and 


^sp,2 


1 

2/<bV^ 


^k^A^Rconv  +  tanh(Atb) 
1  +  A/<bAb^conv  tanh(Atb) 


A  can  be  derived  by 


(14) 

(15) 

(16) 

(17) 

(18) 


7T3/2  1 

V^G 

The  material  thermal  resistance  at  the  base  of  heat  sink  is 


(19) 


Rm  kX  kh(WbLb)  (20) 

When  there  is  a  temperature  difference  between  the  fluid  and 
solid  surface,  according  to  the  Newton’s  Law  of  cooling,  a  convec¬ 
tional  thermal  resistance  needs  to  be  applied.  It  can  be  expressed  as 


^conv 


1 

hAeff 


(21) 


Fig.  3.  One  dimension  thermal  resistance  model  of  the  exhaust  waste  heat  recovery  where  the  heat  transfer  coefficient  h,  according  to  [15],  can  be 
module.  expressed  as 
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Fig.  4.  One  dimension  thermal  resistance  model  of  the  radiator  waste  heat  recovery 
module. 


h  = 


0.037Re4/5Pr1/3 

Hf 


and  the  effective  area  is 

Aeff  =  VfX  (2NWHf)  +  (N-  1  )bW 


where 


Vf  = 


tanh(mHf) 


mHf 


m  = 


h  X 

2 

x  I 

(,,+i) 

I 

l<r  1 

(« 

X  Lj 

I 

2  h 
kf-if 


Thermal  resistance  of  TEG  itself  can  be  expressed  as 
Th-Tl 


rteg  = 


Qh  -  Qc 


(22) 


(23) 


(24) 


(25) 


(26) 


0)  1.6  — 


c 

0) 


fl 
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Fig.  6.  The  performance  of  H-type  TEG  at  circumstance  of  TH  m  340  °C  and 
Tl  -  50  °C. 


2.4.2.  Thermal  resistance  model  of  the  waste  heat  recovered 
from  a  radiator 

Fig.  4  illustrates  a  thermal  resistance  model  using  TEG  to  recover 
waste  heat  from  a  radiator.  This  model  consists  of  the  hot  side,  the 
cold  side  and  TEGs  itself.  The  surrounding  fluid  in  the  hot  side  is  the 
coolant  (water)  of  radiator,  and  air  is  in  the  cold  side.  The  coolant 
temperature  increases  with  the  engine  speed.  A  micro  channel  heat 
sink  is  used  to  simulate  the  function  of  a  radiator.  The  thermal 
resistance  of  the  hot  side  is  composed  of  interface  thermal  resis¬ 
tance,  constriction  thermal  resistance,  convectional  thermal  resis¬ 
tance  and  flow-rate  thermal  resistance.  The  thermal  resistance  of 
the  cold  side  is  the  same  with  the  previous  model. 


100 


80 


ro 

8 

X  40 


20 


0 


Fig.  5.  The  amount  and  percentage  of  heat  loss  from  the  heater  at  different  surface 

temperature.  Fig.  7.  The  variation  of  output  voltage  with  temperature  difference  on  TEG. 
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Fig.  8.  The  performance  curves  of  T-type  TEG  at  circumstance  of  TH  =  150  °C  and 
Tl  =  30  °C. 


When  the  heat  energy  moves  from  a  bigger  area  to  a  smaller 
area,  there  is  a  constriction  thermal  resistance  [16].  It  can  be 
expressed  as 


ftzons 


st  +  O.50r(l  -  e)3/20 
kbira 


where 


(27) 


Table  3 

Cost  evaluations  of  T-type  TEG. 


TGM-127  output  power 

TH  =  90°  C, 

Tl  =  30°C 

Th  =  120°C, 

Tl  =  30°C 

Th  =  150°C, 

Tl  =  30°C 

Maximum  power  (W) 

0.05 

0.08 

0.12 

Power  per  unit  area  (W/m2) 

31.3 

50.0 

75.0 

Cost  per  unit  power  (USD/W) 

920.6 

575.3 

383.6 

When  the  cooling  water  flows  through  the  micro  channel  heat 
sink,  the  convectional  thermal  resistance  at  the  interface  of  water 
and  walls  of  heat  sink  is 


ft 


w,conv  —  t-T  — 


hA  heq(WbbLbb) 


(35) 


According  to  [17],  when  the  heat  flow  from  the  base  of  heat  sink 
to  the  bulk  fluid,  the  equivalent  heat  transfer  coefficient  can  be 
expressed  as 


heq  — 


fim[2JVi,fHf  +  Whib-(N 


—  (N  —  l)fi 


Whib 


(36) 


hm 


Num/<f 

Dh 


(37) 


,  2  WchHf 

h  =  7 - v 

(Wcb+Hf) 


(38) 


The  mean  Nusselt  number,  by  the  discussion  of  Copeland  [18], 
can  be  expressed  as 


tanh(Ar)  +  A/Bi 
1  +  ^/Bi)  tanh(Ar) 


(28) 


Num 


2.22 


RePrDb\°33 
,  Lh,b  ) 


(39) 


A  = 


1 

£y/7T 


(29)  n  =  (8.31G  —  0.02)3 


(40) 


Bi  = 


1 

nkbb(Rconv  + 


r  _  Ar2  +  1 
(30)  _  (Ar  +  l)2 


(41) 


a  = 


b  = 


WhbL} 


h,bLh,b 

7T 


(31) 

(32) 


Ar  =  (42) 

Hf 

When  the  coolant  flows  through  the  heat  sink  channels,  its  inlet 
temperature  will  be  different  with  the  outlet  temperature,  hence 
the  flow-rate  thermal  resistance  needs  to  be  considered.  The  heat 
transfer  rate  is 


£ 


a 

b 


WhbL 


h,bMi,b 


(33) 


T 


b 


(34) 


Table  2 

Cost  evaluations  of  H-type  TEG. 


HZ-2  output  power 

Th  =  300°C, 

Tl  =  10°C 

Th  =320°C, 

Tl  =  30°C 

TH  =  340°  c, 

Tl  =  50°C 

Maximum  power  (W) 

0.35 

0.40 

0.43 

Power  per  unit  area  (W/m2) 

416.2 

475.6 

511.3 

Cost  per  unit  power  (USD/W) 

173.8 

152.1 

141.5 

Table  4 

Performance  coefficients  a,  Re  and  K  of  the  H-type  TEG. 


HZ-2  Performance  coefficient 

a(V/°C) 

Rdn) 

/C(W/K) 

Th  =300°C  Tl  =  10°C 

0.007 

0.05A 

3.77 

0.097 

0.35A 

3.23 

0.70A 

2.92 

Th=320°C  Tl  =  30°C 

0.008 

0.05A 

5.05 

0.097 

0.35A 

3.71 

0.70A 

3.36. 

Th  =340°C  Tl  =  50°C 

0.009 

0.05A 

5.33 

0.096 

0.35A 

4.09 

0.70A 

3.60 

Average 

0.008 

3.89 

0.097 
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Table  5 

Performance  coefficients  a,  RG  and  I<  of  the  T-type  TEG. 


TGM-127  Performance  coefficient 

«(V/°  C) 

RdQ) 

/C(W/K) 

□- 

Th  =  90°C  Tl  =  30°C 

0.011 

0.05A 

2.00 

0.464 

1.24 

—A 

0.35A 

1.88 

0.70A 

2.24 

Th  =  120°C  TL  =  30°C 

0.010 

0.05A 

2.37 

0.309 

1.23 

_ 

0.35A 

2.62 

0.70A 

3.36. 

1 

- 

Th  =  150°C  Tl  =  30°C 

0.010 

0.05A 

3.75 

0.231 

k. 

0 

5  1.22 

_ 

0.35A 

2.90 

o 

0.70A 

2.97 

- 

Average 

0.010 

2.62 

0.334 

1.21 

Q.  —  m*Cp(T0  —  Tjn)  — 


1/m-Cp 


(43) 


In  order  to  define  the  properties  of  coolant  inside  the  heat  sink, 
mean  temperature  is  adopted 


Kfi  = 


_  ^in  +  To 
ean  2 

The  flow-rate  thermal  resistance  is 
1  1 


2-m-Cp  2.(p.Vf).Cp 


(44) 


(45) 


2.4.3.  Operation  parameters  of  TEG 

In  order  to  predict  the  output  behaviors  of  a  TEG,  the  overall 
thermal  resistance  and  the  overall  heat  transfer  coefficient  need  to 
be  known. 
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Fig.  10.  The  variation  of  simulated  output  power  and  thermal  efficiency  with  engine 
speed. 


Qh  -  UhAh(Ts  -  Th) 
Qc  =  UlA l(Tl  -  TC) 


(48) 

(49) 


Equate  equation  (1)  with  equation  (48)  and  equation  (2)  with 
equation  (49),  the  hot  junction  temperature  (TH)  and  the  cold 
junction  temperature  (7L)  can  be  obtained  by 


Uh  =  n - L^~ 

^overall, H^H 

(46) 

Uc  =  R  ~  A 

^overall,  L^L 

(47) 

The  rates  of  supply  heat  (Qh)  and  removal  heat  (Qc)  can  be 
expressed  as 


UhAhTs  +  I<Tl  +  \l2Rc 

(50) 

aI  +  K+  UhAh 

UcAlTc  +  I<Th  +  1/2Rc 

(51) 

UcAl  +  K- Oil 

Fig.  11.  The  simulated  results  of  two  temperatures  (TH,  TL)  related  to  the  coolant 
temperature. 
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Fig.  12.  Performance  of  an  H-type  TEG  at  different  coolant  temperature. 

The  values  of  TH  and  7L  are  derived  by  the  iterative  method. 
Substitute  TH,  TL,  a  and  Rq  into  equation  (11)  to  derive  J0.  Replace 
Th,  Tl>  ol,  Rq,  k  and  J0  into  equation  (1 ),  (2),  (7)  and  (8)  to  find  out  the 
values  of  Qh,  Qo  Rout  and  VW. 

3.  Results  and  discussion 

3.1.  Analysis  of  the  heat  loss 

To  analyze  the  performance  of  a  TEG,  it  is  necessary  to  know  the 
amount  of  heat  delivered  from  the  heater.  Since  the  produced 
power  is  known  from  the  specification  provided  by  the  vendor,  the 
energy  loss  during  the  delivery  from  heater  to  TEG  has  to  be 
measured  first.  Fig.  5  shows  the  amount  of  heat  loss  and  the  heat 
loss  percentage  at  different  heating  surface  temperatures.  Experi¬ 
ment  results  indicate  that  the  heat  loss  increases  with  the 
temperature  of  heating  surface,  but  the  percentage  of  heat  loss 
reacts  in  the  opposite  way.  The  maximum  heat  loss  is  66.7  W  at 
230  °C  and  the  minimum  is  38.6  W  at  50  °C.  The  maximum  heat 
loss  percentage  is  16.3%  and  the  minimum  is  4.9%.  To  have  a  more 
precise  prediction,  the  heat  loss  is  necessary  to  be  considered  in  the 
model. 

3.2.  Performance  curve  and  coefficients  of  a  TEG 

Two  different  TEGs,  H-type  and  T-type,  were  tested  as  baseline 
data.  Fig.  6  presents  the  variation  of  performance  for  an  H-type  TEG 
with  current.  The  output  voltage  decreases  with  the  increase  of 
output  current.  At  the  current  of  0.35  A,  H-type  TEG  produces  the 
maximum  power  of  0.43  W  and  reaches  the  highest  thermal  effi¬ 
ciency  of  1.48%.  Fig.  7  shows  the  relationship  between  voltage  and 
temperature  difference  on  both  sides.  The  voltage  variations  are 
small  at  any  temperature  between  3  different  test  samples. 

Fig.  8  presents  the  test  results  of  a  T-type  TEG  related  to  output 
current,  the  poor  performance  can  be  ascribed  to  the  working 
environment.  The  required  hot  side  temperature  of  H-type  is  higher 
than  the  T-type.  Besides,  H-type  TEG  can  tolerate  larger  tempera¬ 
ture  difference  and  hence  generate  a  higher  output  power.  Tables  2 
and  3  show  the  cost  evaluations  of  two  TEGs.  The  costs  of  TEGs  were 
adopted  with  local  dealers’  prices  in  Taiwan.  Tables  4  and  5  present 


the  values  of  performance  coefficient  a,  R  and  k  under  different 
working  circumstances.  Since  the  performance  disparity,  H-type 
TEGs  were  chosen  for  further  experiments  hereafter. 

3.3.  Modeling  results  of  the  module  recovered  waste  heat 
from  exhaust  pipe 

Fig.  9  illustrates  the  calculated  results  of  TH  and  TL  at  different 
engine  speeds.  The  temperatures  of  TEG  at  hot  and  cold  junctions 
both  increase  with  the  increasing  of  engine  speed.  7H  depends  on 
the  temperature  of  exhaust  gas,  hence  it  is  also  related  to  the  engine 
speed.  The  forced  convection  on  the  cold  side  will  be  improved  by 
increasing  the  engine  speed.  However,  the  thickness  of  TEG  is  only 
5.08  mm,  when  the  heat  spread  from  the  hot  junction  to  the  cold 
junction,  the  cold  side  temperature  is  affected  by  the  effect  of  heat 
conduction  from  the  hot  side.  Therefore,  the  temperature  incre¬ 
ment  in  the  hot  junction  will  be  higher  than  in  the  cold  junction, 
the  temperature  difference  of  TEG  still  increases  with  the  engine 
speed.  The  maximum  temperature  difference  from  simulation  is 
127.4  °C  at  3600  rpm  engine  speed. 

Fig.  10  shows  the  performance  of  a  TEG  at  different  engine 
speeds.  Because  the  temperature  difference  increases  with  the 
engine  speed,  both  output  voltage  and  current  become  higher. 
Therefore,  the  output  power  and  thermal  efficiency  of  TEG  are 
improved  too.  According  to  the  mathematic  model,  the  maximum 
values  of  output  power  and  thermal  efficiency  are  1.24  W  and 
0.396%  respectively  at  3600  rpm  engine  speed. 

3.4.  Modeling  results  of  the  module  recovering  waste  heat 
from  radiator 

Fig.  11  illustrates  the  calculated  results  of  TH  and  7L  of  a  TEG  at 
different  coolant  temperatures.  The  results  show  that,  the 
temperatures  on  both  sides  of  TEG  increase  with  the  coolant 
temperature.  From  the  model,  the  maximum  temperature  differ¬ 
ence  is  42.3  °C  at  the  coolant  temperature  of  90  °C. 

Fig.  12  shows  the  calculated  output  power  and  thermal  effi¬ 
ciency  of  a  TEG  at  different  coolant  temperatures.  The  output  power 
increases  as  the  coolant  temperature  increases.  The  maximum 
output  power  and  thermal  efficiency  is  0.135  W  and  0.135% 
respectively  at  coolant  temperature  of  90  °C. 

4.  Conclusion 

This  research  has  established  fundamental  analyses  of  TE 
modules  applied  on  an  automobile.  The  one  TEG  module  can  be 
extended  to  a  TEG  array  to  produce  more  electricity.  It  is  practical  to 
apply  TEGs  on  an  automobile  to  recover  waste  heat.  The  simulation 
results  of  mathematic  model  have  been  verified  with  experimental 
data  and  showing  the  consistence.  The  characteristic  coefficients  of 
a  TEG  were  derived  from  a  basic  TEG  performance  test.  Under 
designed  working  condition,  a  maximum  power  of  0.43  W  was 
generated  at  0.35  A  current,  and  the  maximum  power  density  was 
51.13  mWcrrr2. 

In  the  simulation  models  of  applications  on  an  automobile,  the 
temperature  difference  between  the  hot  and  cold  junctions  of  TEG 
increased  as  the  engine  speed  or  the  coolant  temperature  increase. 
The  output  voltage,  according  to  the  Seebeck  effect,  also  increased 
as  the  temperature  difference  increase.  Therefore,  the  output 
power  and  thermal  efficiency  can  be  improved.  The  temperature  of 
cooling  water  is  much  lower  than  the  exhaust  gas,  which  caused 
a  smaller  hot  temperature  on  the  TEG  hence  induced  a  poor 
performance.  Different  materials  suitable  for  lower  temperature 
operation  are  needed  for  the  radiator  module. 
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Nomenclature 


A 

area  (m2) 

0 

b 

tunnel  width  (m) 

o,b 

C 

cooling  source 

0,S 

Cp 

specific  heat  at  constant  Pressure  (JKg-1-I<-1) 

P 

D 

diameter  (m) 

S 

H 

high  (m) 

S 

h 

heat  transfer  coefficient  (Wm-2oC-1) 

sp 

I 

output  current  (A) 

T 

I< 

thermal  conductance  (WI<-1) 

W,  ( 

k 

thermal  conductivity  (Wm-loC-1) 

L 

length  (m) 

N 

tunnel  number 

Ref 

P 

output  power  (W) 

V 

perimeter  (m) 

[i] 

Q 

heat  transfer  rate  (W) 

Qh 

hot  side  heat  transfer  rate  (W) 

Qc 

cold  side  heat  transfer  rate  (W) 

[2] 

Rg 

internal  electrical  resistance  of  TEG  (Q) 

R 

thermal  resistance  (KW-1) 

[3] 

S 

heat  source 

T 

temperature  (°C) 

[4] 

t 

thickness  (m) 

U 

overall  heat  transfer  coefficient  (Wm-2oC-1) 

V 

voltage  (V) 

[5] 

W 

width  (m) 

[6] 

Dimensionless  number 

[7] 

Bi 

Biot  number  Bi  =  .  .  ,n - =— 

7r/<bb(RConv  +  Rfi) 

Nu 

Nusselt  number  Nu  =  /i*L//<fluid 

[8] 

Pr 

Prandtl  number  Pr  =  p»Cp/k 

Re 

Reynolds  number  Re  =  VHf  j v 

[9] 

Greek 

symbols 

a 

Seebeck  coefficient  (V°C-1) 

[10] 

P 

dynamic  viscosity  (N*5*m-2) 

V 

kinematic  viscosity  (m2*s-1) 

[11] 

V 

efficiency  (%) 

Subscripts 

[12] 

b,con  base  mean  temperature 

b 

base  of  Teflon  block 

[13] 

C 

cooling  fluid  temperature 

[14] 

c 

contact  surface 

ch 

channel 

[15] 

conv 

convection 

cons 

constriction 

[16] 

eq 

equivalent 

eff 

effective  area 

[17] 

f 

fin 

[18] 

fl 

flow  rate 

G  TEG 

h  hydraulic 

h, b  heat  sink  base 

H  hot  side 

i  interface 

i, b  inner  base  of  Teflon  block 

i,s  inner  side  of  Teflon  block 

L  cold  side 

1  heat  loss 

m  material 

output 

outer  base  of  Teflon  block 
outer  side  of  Teflon  block 
constant  pressure 
heat  source  temperature 
side  of  Teflon  block 
spreading 
Teflon  block 
w,  cv  convection  between  water  and  wall 
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